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Abstract 

We perform a detailed analysis on Affleck-Dine leptogenesis taking into account 
the thermal effects on the dynamics of the flat direction field (p. We find that an 
extremely small mass for the lightest neutrino z^i, niui ^ 10~^ eV, is required to 
produce enough lepton-number asymmetry to explain the baryon asymmetry in the 
present universe. We impose here the reheating temperature after inflation Tr to 
be Tr ^ 10^ GeV to solve the cosmological gravitino problem. The required value of 
neutrino mass seems to be very unlikely the case since the recent Superkamiokande 
experiments suggest the masses of heavier two neutrinos 1^2 and ^3 to be in a 
range of 10~^-10~^ eV. We also propose a model to avoid this difflculty based on 
the Peccei-Quinn symmetry, where the required neutrino mass can be as large as 
niuA — 10"'* eV. 



1 Introduction 



There have been growing interest in leptogenesis to account for the baryon asymmetry 
in the present universe since the Superkamiokande collaboration presented convincing 
evidence of the atmospheric neutrino oscillation P). In fact, various mechanisms p, ^, Q 
for the leptogenesis have been proposed so far. Among them the leptogenesis by fiat 
directions in super symmetric (SUSY) theories, i.e., Affleck-Dine (AD) mechanism 0, 
is expected to work even with relatively lower reheating temperatures Tr (e.g., Tr ^ 10^ 
GeV) and hence it seems cosmologically safer than others, since the number density of 
gravitinos produced in the reheating processes is proportional to the reheating tempera- 
ture Tr of inflation |^.[| However, it has been pointed out p, |10] that the flat direction 
flelds (p start their coherent oscillations at an earlier time than usually estimated due 
to interactions with thermal plasma of inflaton-decay products and that the resultant 
baryon asymmetry is strongly suppressed. 

In this paper we perform a detailed numerical analysis on the AD leptogenesis taking 
into account the thermal plasma effects. It has been argued in Ref. that a relatively 
low reheating temperature Tr < 10' GeV is necessary to escape from the thermal effects. 
We conflrm that the baryon asymmetry is indeed suppressed for Tr ^ 10^-10^ GeV due to 
the thermal effects, when the lightest neutrino mass is m,yi ^ 10^^ eV. On the other hand, 
we flnd that the thermal effects become signiflcant for a higher reheating temperature, 
Tr>W^-10' GeV x(m^i/10-9 eV)-\ if m^i < IQ-^ eV. 

Furthermore, we estimate how much the resultant baryon asymmetry is suppressed 
for such a high reheating temperature by analytical and numerical calculations. What 
we have found is that an ultralight neutrino of mass m^^i ^ 10~* eV is required to obtain 
the desired baryon asymmetry in the present universe, ra^/s ~ (0.1-1) x 10~^°, when we 
impose the reheating temperature Tr to be Tr ^10^ GeV to avoid the the overproduction 
of gravitinos 0. Here, ub and s are the present baryon- number and entropy densities, 
respectively. 

We note, however, that the above neutrino mass is not necessarily the same as the 
mass to be observed today, if effective Majorana masses for right-handed neutrinos in 
the early universe are different from the values in the true vacuum. To demonstrate this 
point we show an explicit model based on the Peccei-Quinn symmetry [|lT|, which is one 
of the most attractive solutions to the strong CP problem. In this model we flnd that 
the mass for the lightest neutrino in the true vacuum can be as large as 10~^ eV. This 



^ The gravitinos may be also produced nonthermally We neglect such an effect in this paper, 
since it depends on the details of inflation models. 
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is very close to the masses of heavier two neutrinos, mj,2,3 — 10 ^-10 ^ eV, suggested 
from the recent Superkaimokande experiments [@, |l^ . 



2 Flat directions for the Leptogenesis 

We adopt fiat directions = Li = (jf / \pl (a family index 2 = 1,2, 3) originally proposed 



by Murayama and one of the authors (T.Y.) |T3[ as the Affieck-Dine fiat direction field 
for leptogenesis, where Li are scalar components of the chiral multiplets Li of SU{2)l- 
doublet leptons. We have effective dimension-five operators in the superpotential, 

W=^ {L.H^f = ^ {L^H^f . (1) 



Here, = {Hy) = sin /? x 174 GeV and we have used the seesaw formula ||T4 



= -M-- (2) 

Hereafter, we take sin/5 ~ 1, and we suppress the family index i for simplicity. As we 
will see later, the relevant flat direction for the most effective leptogenesis corresponds 
to the first family (i.e., (p/\/2 = Li = Hu)- The potential of the flat direction is given 
by 

where m and A are SUSY-breaking mass parameters and the last term is derived directly 
from the superpotential in Eq. ([^). We choose m ~ |A| ~ 1 TeV. 

During inflation the energy density of the universe is dominated by the inflaton x- 
After the end of inflation the inflaton x starts its coherent oscillation and the energy 
density of the universe is dominated by the inflaton x also in this period. The non- 
zero vacuum energy of x induces additional SUSY breaking effects. Therefore, the flat 
direction receives an additional SUSY breaking mass of the order of H in the X" 
dominated era. Here H is the Hubble parameter of the expanding universe. There is 
also an A-term proportional to the Hubble parameter, V ~ [H/M)(j)'^. Thus, we have 
the additional potential is given by 

5V = -chH'IM' + ^ {aH<P' + a*^<P*') , (4) 

= {Hu) / (Hd), where Hu and Hd are Higgs supermultiplets which couple to "up" -type and 
"down"-type quarks, respectively. 
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where ch and an are real and complex constants respectively, and we take ch > 0. (See 
the next section.) Hereafter, we assume ch — \aH\ — 1-Q 

The crucial point in the following discussion is that the decay of the inflaton x occurs 
during the coherent oscillation of while it completes much later than the beginning 
of the oscillation. Thus, even in the x oscillation period, there is a dilute plasma with a 
temperature T = {T'^M^HY^^ {M^ = 2.4 x 10^^ GeV is the reduced Planck mass) arising 
from the inflaton decay , although most of the energy density of the universe is carried 
by the coherent oscillation. Notice that in this dilute plasma the temperature decreases 
as T oc H^^^ rather than as T oc if^/^ in the usual radiation dominated universe. 

It is quite plausible that the fields ipk which couple to the are produced by the 
inflaton decay and/or by thermal scatterings if their effective masses fk\4'\ are smaller 
than the temperature T, 

M0|<T, (5) 

and hence they must be included in the plasma. Here, fk correspond to their Yukawa or 
gauge coupling constants of the fiat direction (we take fk real and positive). 

As stressed in Refs. |T0| the thermal effects from the dilute plasma may affect 
the dynamics of the (p field. In fact, the fiat direction (p receives a thermal mass of 
order ~ f^T if the condition Eq. (|^) is satisfied. Then the total effective potential for 
including the Hubble-induced terms and the thermal mass term is given by |ffOl 



\ fk\<t>\<T 

where Ck are real positive constants of order unity and a = A/m is a complex constant 
of order unity. The summation of fk means that only the fields in thermal plasma [i.e., 
their couplings fk satisfy the condition Eq. (^)] can induce the thermal mass for 0. 

We include all the couplings relevant for the fiat direction Hu = Li = 0/-\/2 in the 
SUSY standard model, i.e., the gauge couplings for the SU{2)l and f/(l)y gauge groups, 
and Yukawa couplings for up-type quarks and charged leptons. The couplings fk, which 
are redefined so that the masses for the fields ipk couple to become /fc|0|, are listed in 
Table.[l| with the coefficients Ck of the thermal mass for 0. Among the gauge supermulti- 

The SUSY breaking effects due to the finite potential energy of the inflaton x depend on the details 
of the Kahler potential. In fact, the coefficients of the Hubble-induced terms, ch and an in Eqs. (^) and 

can be much smaller than order unity. However, we simply assume here ch — \o-h\ — 1 since the 
conclusion does not depend much on these parameters unless ch ^ —1- (See, for example, Ref. pq].) 
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fk (coupling for 0) 


V 




g2 


{a = t,c,u) 
^ 3 




Ck 


i 

4 
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Table 1: The couplings of the flat direction field to other fields and the coefficients Ck 
of the thermal mass of induced by these fields, gi and g2 denote the gauge couplings 
for the U{1)y and SU{2)l gauge group respectively, and i/a (a = t,c,u) and yii are the 
Yukawa couplings for up-type quarks and charged lepton, respectively. 

plets, the fiat direction does not couple to one U{1) group which remains unbroken by 
the condensate of 0. Thus, there are one Z-like and two Pi^-like massive gauge multiplets 
as in the SUSY standard model. For the couplings of the charged leptons, we should 
be careful. The fiat direction receives its thermal mass from only one linear combina- 
tion of charged leptons. The effective Yukawa coupling yi-^ is determined by the Yukawa 
couplings for the charged leptons and the mixing matrix among gauge eigenstates of the 
neutrinos. (Note that the relevant is the fiat direction which corresponds to the first 
family Ui of mass eigenstates.) Taking large mixing angles both for atmospheric and 
solar neutrino oscillations indicated by the recent Superkamiokande data p|, we find 
yi^ ~ (9(0.1-1) X y^. Thus, fk are roughly 10~^-1 in the SUSY standard model. We will 
see in Sec. § that the Yukawa coupling for the up quark induces relevant thermal effects, 
while other couplings turn out to be too large to satisfy the condition Eq. in most of 
the cases. 



3 Evolution of the flat direction 

In this section we discuss the evolution of the flat direction 0, which is described by the 
equation of motion with the potential Kotai in Eq- (H), 

+ 3i70+^ = O, (7) 

where the dot denotes the derivative with time. 

During inflation the Hubble parameter H^^f is much larger than the soft SUSY break- 
ing mass, i^inf ^ m. We assume that the additional mass squared is negative {ch > 0) 
and causes an instability of at the origin, that is, miff = w? — cu^HmiY < 0- Then, 
the minima of the potential are given by 

101 ^ ^MH,^f, (8) 
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- arg(a//) + (2n + l)7r 
arg(0) ^ ^ — — n = 0---3. (9) 

Note that there is no plasma and hence no thermal mass term during the inflation. 
Because curvatures of the potential around the minimum along both the radius and 
phase directions are of the order of the Hubble parameter -ffjnf, the fiat direction runs 
to one of the four minima {n = uq) from any given initial value and is settled down it. 

After the inflation ends, the infiaton x starts to oscillate, and there appears a dilute 
plasma with a temperature T = {T'^M^.HY^^ . Then, the dynamics of the fiat direction 
is determined by the total potential Kotai in Eq. (H). 

In the period of the infiaton oscillation, the minimum of the potential moves toward 
the origin = as the Hubble parameter H decreases. As long as the potential for is 
dominated by Hubble-induced (mass and A) terms and also |0|^ term, the fiat direction 
field always trades the following instantaneous minimum of the potential |P, [15|, 



101 ~ ^MH, (10) 

- arg(a/^) + (2^0 + l)7r 
arg(0) ~ . (11) 

This is because the curvatures around the minimum along both radius and phase direc- 
tions are of the order of H also in this period, and hence the always catches up the 
instantaneous minimum. 

As the universe evolves, the field value |0| decreases with time. When the field value 
|0| satisfies the condition Eq. (|]), the fields ipk which couple to the can enter in the 
dilute plasma so that they induce a thermal mass Cfc/|T^ for 0. If this thermal mass 
exceeds the Hubble-induced mass, i.e., 

H'< ^^flT\ (12) 

h\4>\<T 

the dynamics of the fields is drastically changed. The equation of motion for at this 
epoch is approximately given by 

+ 3i/0 + CkflT^ct) = 0. (13) 

This equation can be solved analytically and the solution is given by 



0(t) = 










01 


[^2/3 ^^Vsl^)] 


+ 02 


[^2/3 "^-^/Sl^)] 








\ 1/4 




z = 


u 









(14) 
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where 0i and 02 are constants and is the Bessel function. We see that the oscillates 
around the origin = 0. The time scale of the oscillation is ~ (/^M^,T|) and the 
amplitude of the oscillation is damped as |0| ~ t~'^'^ ~ 

Now let us estimate the cosmic time t^^^ when the flat direction starts its oscillation 
due to the thermal mass term, since the time t^^^ plays a crucial role in the production 
of the lepton-number asymmetry. First of all, we consider a simple case in which there 
is only one field ipi with a coupling fi. It is easily found from Eqs. and ( |12|) that the 



following two conditions must be satisfied in order to cause the 0's oscillation. 



1 M.TE 

1/3 



H < (clftM^Tl)'\ (16) 



where we have used the relations |0| ^ V MH and T = {T^M^Hy^"^. If the coupling /i 
is small enough, the ipi can easily enter in the thermal plasma, but the induced thermal 
mass for is smaller than the Hubble-induced mass H at the beginning. However, since 
the temperature decreases more slowly than the Hubble parameter H {T oc H^^^) as 
the universe expands, the thermal mass eventually exceeds H . Then the 0's oscillation 
starts at the time when the condition Eq. ( [T^ ) is satisfied. On the other hand, if /i is 
very large, the would-be thermal mass can be large enough to exceed the Hubble mass. 
However, the thermal mass does not appear until the |0| becomes small enough to satisfy 
Eq. (^. In this case the starts to oscillate at the time when the condition Eq. (|I3]) is 
satisfied. Therefore, we find the Hubble parameter at t = t'^^^^ to be 



H± = min 



osc 



IMA , (.;/jM.rlV" 



(17) 



ft M2 

Notice that in order to have a significant thermal effect, as discussed in Ref. |]I0[, a 
coupling fi with intermediate strength is necessary. This is because too large f\ can not 
satisfy the condition Eq. (^), while too small /i can not give a large thermal mass for 
which satisfies the condition Eq. (p!2D. 

It is easy to apply the above discussion to the case of more than one couplings. If 
there is another coupling fi which can satisfy the both conditions Eqs. (|15D and (0) 
earlier than /i, the flat direction starts its oscillation earlier. Therefore, the Hubble 
parameter is given byQ 

1 M,T| / o,4,^^2\V3^ 



^oL = max 



mm 



4 T\/r2 ' y^iJi 



(18) 



We assume here hierarchical couphngs and neglect effects of the summation. 
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Thus, becomes larger (i.e., the oscillation begins earlier) as we take the reheating 
temperature Tr higher. One can see from Eq. ([18| ) that it is nontrivial which cou- 
pling among /j makes the 0's oscillation earliest, i.e., it is not always the coupling with 
strongest or weakest strength. If one takes the reheating temperature higher, a larger 
coupling makes the oscillation earlier and determines the time H^^. 

As the universe becomes cooler, the thermal mass for the becomes eventually 
smaller than the original mass m. Then, the evolution of (p is described by the equation 
of motion 



The flat direction oscillates around the origin with a time scale ~ m~^, and the amplitude 
of the oscillation is damped as \(f)\ ~ t^^ ~ as the universe evolves. 

We have, so far, assumed that the thermal mass term dominates the potential when 
the oscillation of (p starts. However, the expression of in Eq. ( p^Sf ) could be smaller 
than the original mass m, depending on the parameters M and Tr, which conflicts with 
the above assumption. If the thermal mass does not exceed the Hubble mass H until 
H "r^i m, the original mass term m^|(/)p dominates the potential for H^m. Then, the 
thermal effects induced by ipk in the dilute plasma do not affect the evolution of 0. In 
this case the (p field begins its oscillation at if ~ m. 

In summary, the oscillation of starts when the Hubble parameter of the universe 
becomes H ~ iiosc, where 

1. Hose = -f^osc if thermal mass term dominates the potential before H c:^ m, 

2. Hose = rn, if otherwise. 

We show a schematic behavior of i^osc in Fig. |l] for illustration. One can see from 
Eq. (|T8|) that the Hubble parameter i^osc depends on the reheating temperature non- 
trivially. We consider three couplings /i, /n and /m assuming a hierarchy fi <^ /n ^ /m. 
The region of Tr can be divided into three regions I, II and HI as in Fig. [l|, where H^^^ 
is determined by /i, fu and /m, respectively. [See the discussion below Eq. (|18|).] The 
dependence of H^^ on Tr can be understood from Eq. (^) or (|TB]). For example, if Tr 
lies in the region II, Hose behaves as 



+ 3H^ + = 0. 



(19) 




for Tr < Ten 
for Tr > Ten 



(20) 



where the critical temperature Ten is 




(21) 
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Figure 1: A schematic behavior of the Hubble parameter ifosc when starts its oscilla- 
tion. The oscillation of starts at H^sc = in region la, while the oscillation starts at 
an earlier time (-ffosc = -f^osc) regions lb - III. 



We can also see the same behavior in the regions I and III. Furthermore, when the 
reheating temperature Tr is low enough (i.e., in the region la), the thermal effects 
become negligible and the oscillation of begins at i^osc = We see from Fig. that 
the oscillation time -ffosc, which is crucial for the estimation of the lepton asymmetry 
(see the next section), depends on the reheating temperature and also the spectrum of 
the couplings in a very complicated way. 



4 Lepton asymmetry 

Now we are at the point to calculate the lepton asymmetry produced by the field. 
Since the field carries the lepton charge, its number density is related to the lepton 
number density as 
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First we consider the case ifosc = -f^osc > Suppose that the original A-term would 
vanish, i.e., A = ma = 0. If this is the case, the flat direction is always trapped in 
one of the valleys induced by the Hubble A-term during the both periods t < t^^^ and 
t > t*gj,, and the direction of the valley does not change with time. Therefore, there is 
no force which causes the motion of (p along the phase direction and no lepton-number 
asymmetry is produced. However, we have the original A-term, and the phase of (j) is 
kicked by the relative phase difference between a and an in Eq. (^). The phase of 
changes during its rolling towards the origin because the Hubble parameter H decreases, 
and the direction of the true valleys changes with time. Therefore, the original A-term 
[the first term in Eq. plays a role of the source of the lepton asymmetry. 

One might wonder if the Hubble A-term [the second term in Eq. (P5D] gives a larger 
contribution to the lepton asymmetry since if ^ m in the present situation. However, 
since the (p almost traces one of the valleys determined mainly by the Hubble A-term, 
Im(aH0^) in Eq. ( p3D is highly suppressed compared with Im(a0''). In fact, we have found 
numerically that the contribution of the second term in Eq. (^) is always comparable 
or less than that from the original A-term. Thus, we neglect the lepton asymmetry 
produced from the second term in Eq. (|23| ) in our analytic calculation, for simplicity. By 
integrating Eq. we obtain the resultant lepton number at the time t as 



where R denotes the scale factor of the expanding universe, which scales as i? oc H^"^^^ oc 
t^/^ in the universe dominated by the oscillating infiaton x- Notice that the integrand of 
the first term in Eq. (pSD is almost constant, since 101 ~ y/MH and R oc R-"^/^. On the 
other hand, the second term gives only small contribution to the total lepton asymmetry 
because of the following two reasons; (i) Im(a0^) changes its sign rapidly due to the 
oscillation of 0, and (ii) even if the sign would remain unchanged, the amplitude of 0^ 
is damped as |0|^ ~ (see Sec. ||). 

Therefore, the production of lepton asymmetry is suppressed for t ^ t^^^, and the 




(24) 



(25) 
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resultant lepton asymmetry at t ~ t]^^ is given approximately by 



.th 

''osc 



= ^mMHli6,s, (26) 

where 6es = sin(4arg0(t*gj.) + arga) represents an effective CP violating phase and we 
assume 6es — Notice that we have used the fact |0(tosc)l — \J ^^^c- 

The lepton-to-entropy ratio when the reheating process of inflation completes (T = 
Tr) is estimated as 



MTji ( m 



s 12M? \ H± 



4ff. (27) 



This ratio takes a constant value as long as an extra entropy production does not take 
place at a later epoch, since both ul and s decrease at the same rate as the universe 
expands. In deriving Eq. (^) we have assumed that the oscillation starts before the 
reheating process of inflation completes. 

Next we turn to the case without thermal effects, i.e., i^osc = '^tt-- Detailed analysis 
for this case is found in Refs. |P, 15|. Even in this case, the production of the lepton 



asymmetry becomes ineffective after the (p starts the oscillation (for H <^ i^osc = "o^)- 
Thus, the lepton asymmetry can be estimated by replacing with m in Eq. (^T]) as 



T - T2Mf^'=^- ^^^^ 
The produced lepton asymmetry is partially converted [ffl] into the baryon asymmetry 



through the "sphaleron" effects |17], since it is produced before the electroweak phase 



transition at T ~ 10^ GeV. Then, the present baryon asymmetry is given by |T8| p| 

^ = A!^. (30) 

s 23 s ^ ^ 



Now we have formulae of the lepton asymmetry Eqs. (|27|) and (|28|) . Notice that an 
additional suppression factor m/H^^ (-^osc > appears in Eq. (|27|) , compared with 

^ In the present analysis we neglect the sign of the produced lepton asymmetry. To be exact, the 
lepton asymmetry should be negative, since the present baryon asymmetry is given by p8{ 

Ub _ 24 + iNn riL 

rather than Eq. (pO|). Here Nh denotes the number of the Higgs doublets. Then, ub/s — — (8/23)nz,/s 
in the minimal SUSY standard model with Nh = 2. 
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Figure 2: A schematic behavior of the produced lepton asymmetry ul/ s. Regions I, II, 
III correspond to those in Fig. |l|. 

Eq. (pHD- This suppression clearly represents the effect of the early oscillation due to 
the thermal effects of the dilute plasma pointed out in Ref. 0, 0. It should be also 
noted that the produced lepton asymmetry becomes larger as the scale M increases, i.e., 
as the neutrino mass m,^ becomes smaller [See Eqs. (0), (0) and (PB])]. Therefore, the 
leptogenesis is most effective for the fiat direction corresponding to the lightest neutrino 
z/i, as mentioned in Sec. ^ 

The produced lepton asymmetry depends on the reheating temperature Tr. We show 
again a schematic behavior of the lepton asymmetry in Fig. |^, where regions I, II, and 
III correspond to those in Fig. [1|. If the reheating temperature Tr is low enough, i.e., 
there is no thermal effect (ifosc = the lepton asymmetry ul/s is just proportional 
to Tr as in Eq. (^, 

- « n\ (31) 

s 

However, it is not the case when the reheating temperature becomes higher (i.e., iTosc = 
since H]^^ depends on Tr (see Fig. 0). Depending on the condition which deter- 
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mines H^^, the lepton asymmetry behaves as 

^ oc T^'orT'J'. (32) 

The exphcit formulae of the produced lepton asymmetry are easily derived from Eqs. (0) 
and (^). For example, if the reheating temperature lies in the region II, the resultant 
lepton asymmetry is given by 



s l2M^Tn 



for Tr < Ten, 

for T^>Tc„, (33) 

where Ten is given by Eq. (pi]). 

Now, we discuss how much lepton asymmetry is generated in the SUSY standard 
model. We have performed a numerical calculation to follow the evolution of the flat 
direction for given M and Tr, including all the terms in the potential in Eq. (^. As 
for the thermal mass term, we have included all the couplings listed in Table.|l|. In Fig. ^ 
we show the produced lepton asymmetry for the lightest neutrino mass Tfiui — 

10-^ eV. 

Here, we take m = 1 TeV, tan/5 = 3 and aig{a/aH) = 7r/3 (we take this parameter set, 
hereafter). As shown by the dotted line in Fig. 0, the produced lepton asymmetry is 
proportional to the reheating temperature Tr, if one neglects the thermal effects (i.e., 
takes all the coupling fk = 0). In this case, we can derive from Eq. (^) the formula for 
the lepton asymmetry 

which is in a very good agreement with the result of the numerical calculation in Fig. ^ 
The above expression Eq. (|33) is also consistent with that obtained in the previous 



works 1^, |T5| . Notice that the produced lepton asymmetry is proportional to the reheating 
temperature and inversely to the lightest neutrino mass. Therefore, in order to explain 
the baryon asymmetry in the present universe a high reheating temperature of Tr ^ 10^ 
GeV is required for mj.^ ~ 10~^ eV. However, for such a high reheating temperature 
the effects of the thermal mass on the dynamics of is significant, and hence the 0's 
oscillation in fact begins at i^osc = H^sc rather than i^osc = Thus, the formula Eq. (0) 
can not be applied to the case of m,yi ~ 10"^ eV for the high reheating temperature. 

The solid line in Fig. ^ represents the lepton asymmetry with the thermal effects. 
Even in this case, the produced lepton asymmetry is linear in Tr for a low reheating 
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Figure 3: The produced lepton asymmetry rii/s for niyi = 10^^ eV. The sohd hne 
represents the lepton asymmetry including the thermal effects. The dotted line represents 
the lepton asymmetry when one neglects the thermal effects (i.e., /fc = 0). We take m = 1 
TeV, tan/5 = 3 and arg(a/aff) = tt/3. 



temperature region of Tij < 4 X 10^ GeV as in Eq. (|3^) . However, we can see that the 
thermal effects become significant and the lepton asymmetry is suppressed for Tr ^ 10^- 
10^ GeV, which is consistent with the argument in Ref. 0. This is because the thermal 
mass term makes the 0's oscillation earlier (-ffosc = ^ m), which reduces the lepton 
asymmetry. In fact, the behavior of the produced lepton asymmetry shown in Fig. ^ can 
be understood by our analytic formula Eq. (|^) and Fig. ^. 

Fig. H shows the contour plot of the produced lepton asymmetry rii/s in the M-T^ 
plane. We show both results by the numerical calculation and also from the analytic 



formulae Eqs. (^) and (|28|) . We confirm that the analytic estimation discussed above 
reproduces very well the result obtained by the numerical calculation. 

In the lower reheating temperature region where ifosc = the lepton asymmetry 
Ul/s is proportional to MTr [see Eq. (|^)]. On the other hand, in the region where 
^osc = ^osc > ^) ^l/s is proportional to M^T^^ or MT^^^ [see Eqs. (|18D, and 
(0)], depending on the condition which determines the cosmic time tosc when the early 
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M [GeV] 



Figure 4: The contour plot of the lepton asymmetry n^/s in the M-Tr plane. The solid 
lines show the contour lines of rii/s obtained by numerical calculation. The dashed lines 
show the contour lines of rii/s obtained from analytic formulae. Corresponding values 
of ul/s are also represented. We take m = 1 TeV, tan/? = 3 and arg^a/an) = 7r/3. 

oscillation of starts. 

1 /3 

In most of the region in Fig. | the lepton asymmetry behaves as n^/s oc MT^ and 
the thermal effects are controlled by only the Yukawa coupling constant for the up quark, 
which is roughly of order 10~^. Only in the parameter region M < lO^i GeV and Tr > 10^ 
GeV, larger couplings (say, the Yukawa coupling for the charm quark) become effective. 
Note that the result shown in Fig. ^ is almost independent of tan j3 as long as sin /5 ~ 1 
since the Yukawa coupling for the up quark is Uu = rriul {H^) = m^J/(174GeV x sin/3), 
where m„ is the mass for the up quark. 

The crucial observation is that the resultant lepton asymmetry is suppressed by the 
thermal effects when the reheating temperature is high enough. The critical value of the 
reheating temperature, Tfjcr, is estimated analytically from Eq. (^) as 




(35) 
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for the region of M < Mcr, and 

Tn.. = /^-^ = 4 X 10= GeV [^^^) (^) . (36) 

for the region of M > Mcr, where M^r is given by 

= -\rm = 4 X 10^2 GeV ( . (37) 

This critical value for the scale M corresponds to the lightest neutrino mass of rrtyi = 
0.9 X 10~^ eV for m = 1 TeV. Here we have used the fact that only the up-quark Yukawa 
coupling is effective for the thermal effects. Therefore, we find that the thermal effects 
from the dilute plasma is suppressed in the region Tr < lO^-lO^ GeV for m^i > IQ-^ eV. 
This is consistent with the result obtained in Ref. 0. However, for the lighter neutrino 
mass region m^i ^ 10~^ eV, this upper bound on the reheating temperature in order 
to neglect the thermal effects increases as M (i.e., inversely proportional to the lightest 
neutrino mass m^i). In fact, it helps us a lot with the discussion in the next section. We 
also confirm the result by a numerical calculation. 

In this section, we show that the thermal effects on the AD leptogenesis is very 
significant. In fact, when we require the reheating temperature should be Tr ^ 10^ GeV 
to avoid the gravitino problem 0, we need very a large value of M, say, 3 x 10^^ 
GeV, in order to obtain the desired lepton asymmetry rii/s ~ 10^^'^-10~^. This large 
value M corresponds to an extremely small neutrino mass rriyi < 10"*^ eV. 

5 Large effective masses for right-handed neutrinos 

Let us discuss the origin of the operators in Eq. (|l]) in the presence of the heavy right- 
handed Majorana neutrinos Nj (/ = 1,2,3). The relevant superpotential is 

W = Kj^,NjL,Hu + ^MnjNjNj , (38) 

where we have used a base where the mass matrix for the heavy neutrinos Nj is diagonal. 
Then, through the see-saw mechanism the heavy Majorana neutrinos Nj induce the 
effective operators in Eq. (|l|) at low energies, 

W = l^^W.H. . (39) 



These operators gives light neutrino masses as 

T 



m.k, = vl^^ . (40) 
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It is quite natural to consider that the Majorana masses M^j for Nj are given by a 
vacuum expectation value (vev) of some field X, which is a singlet under the standard- 
model gauge groups, with the superpotential 

W = hiMiXNiNi . (41) 

If it is the case, masses of the heavy Majorana neutrinos Nj in the early universe may be 
different from those in the present universe. Namely, the extremely small neutrino mass 
discussed in the previous section {ni^i ^ 10~^ eV) may not be the same as the observed 
mass. In particular, when the field X is responsible to the Peccei-Quinn symmetry 
breaking fill], we have necessarily a flat direction containing the X field. Thus, the X 
field may have a large value as X ~ M* during the inflation whereas it has a vev of 
(X) ~ Fa in the true vacuum. Here, the Peccei-Quinn breaking scale Fa is constrained 
by laboratory experiments, astrophysics, and cosmology as Fa — 10^''-10^^ GeV ||19|| .p] 



In this case, neutrino masses in the early universe are much smaller than those observed 
today (they are suppressed by the factor X/Fa — M^/Fa — 10^-10*^). Therefore, the 
extremely small neutrino mass required for a successful leptogenesis may be naturally 
explained in this scenario. 

To demonstrate our point, we consider the following superpotential for the Peccei- 
Quinn symmetry breaking sector, 

W = \Y{XX-Fl) , (42) 

where A is a coupling constant, and F, X, and X are supermultiples which are singlets 
under the standard-model gauge groups and have 0, +1, and —1 Peccei-Quinn charges, 
respectively. From Eq. (^21) we see that there exists a fiat direction XX = F^. We 
parameterize this direction by the scalar field a called "saxion". This saxion a receives 
a soft SUSY breaking mass of mo- — WI3/2 (^3/2 is the gravitino mass), and the X and X 
have vevs of the order of Fa if the soft SUSY breaking masses for X and X are almost 
the same. {C = m^|Xp + m^lXp with nix — ""^x — ^3/2)- 

Along this fiat direction the X field may have a large value during the inflation, if the 
additional SUSY breaking effects induce a negative mass squared for X. In this case the 
supergravity effects prevent the X field from running over M^, and give an initial value 
X ~ M^, We here assume X ~ M^, and hence X ~ Fa/M^ at the end of infiation.[| 



^ The upper bound on Fa is relaxed as Fa ^ 10^^ GeV, when there is an extra entropy production 
after the QCD phase transition [ po| . 

In this case, the isocurvature fluctuation associated with the axion is very suppressed . 
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We turn to calculate the produced lepton asymmetry in the presence of the X field. 
We see from Eq. (0) that the masses for Ni are 



Mr, = hNiX . (43) 

Then, we find 



= -^Kfjhiyj ^h^j j . (44) 

The vev of X in the true vacuum, (X) ~ Fa, is supposed to give the neutrino masses 
observed today. As discussed in the previous sections, the lightest neutrino ui {rriyi <^ 
m;^2 ^ "^vi) is relevant for the AD leptogenesis. 

As for the neutrino masses we adopt the recent data from the Superkamoikande ex- 
periments. The atmospheric neutrino oscillation is well explained by the Vpri^r oscillation 
with m,^3 ~ 5 X 10~^ eV 0. On the other hand, the large MSW solution to the solar 
neutrino problem is favored by the recent Superkamiokande data [Q, which indicates 
that m^2 — 5 X 10~^ eV. However, it is very difficult to determine the lightest neutrino 
mass niui- Here, we simply take m^i = 1 x 10~^ eV for a representative value. 

Then, the scale M in the previous sections is estimated as 

Notice that as shown in Fig. ^ such a scale M for (X) ~ Fa can not produce sufficient 
lepton asymmetry to explain the baryon asymmetry in the present universe. However, 
the large initial value of X {Xq ~ M^) gives 



M = 7 X 10^^ GeV 



10-^ eV\ /I0i° GeV\ f Xq 



which is translated into the lightest neutrino mass in the false vacuum 
/ 1-? / ^ul \ ^ Fn \ f 



Thus, such a high scale M in Eq. (^6]) can give a neutrino mass much smaller than that 
observed today. It should be noted that the fiat direction saxion cr stays at Xq ^ until 
H ^ rricr — Tns/2 due to the friction of the expansion of the universe. Furthermore, we 
find from Eqs. ( pS] ) and ( |36D that for such a high scale of M ~ 10^^ GeV we can estimate 
the lepton asymmetry without the thermal effects even if the reheating temperature is 
as high as Tr ~ 10^ GeV. 
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Therefore, for the reheating temperature of ^ 10® GeV, the produced lepton asym- 

as 



metry is estimated by using Eq. 

s 12M2 I m,iF„ 



off 



- ' ^ (lo^vj [Wj ■ (^«) 

This shows that too large amount of the lepton asymmetry is produced through the AD 
mechanism. However, this asymmetry is sufficiently diluted, since there exists substantial 
entropy production at low energies. The saxion a is produced in the coherent oscillation 
after the lepton asymmetry is frozen, and the oscillation energy dominates the universe.0 
Then, the late decay of the saxion increases the entropy of the universe and dilutes the 
lepton (baryon) asymmetry substantially. 

The saxion begins the coherent oscillation at if ~ m^r — with the initial ampli- 
tude Xq ~ M^,. The oscillation energy at that time is given by p^- ~ m^XQ/2. Note that 
the saxion oscillation starts before the reheating process of the inflation takes place.0 
Then, when T = Tr, the ratio of p„ to the entropy density of the universe is estimated 
as 

^ = -Tn(^)\ (49) 

Notice that p„ decreases at the rate R^^ as the universe expands, while the radiation 
energy density decreases as Therefore, the oscillation energy of the saxion domi- 

nates the energy of the universe soon after the reheating process completes. This energy 
is transferred into thermal bath by the saxion decay. The saxion decays into two gluons 
with the partial decay rate 

rv Tn 

Through this decay the universe is reheated again, and its reheating temperature Tg- is 
estimated as 



^ The saxions are also produced by the thermal scatterings in the reheating period. However, the 
energy density of the saxions produced thermally is much smaller than that produced in the coherent 
oscillations, and hence we neglect it. 

9 It is the case for Tr<2x 10^" GeV(TO^/l TeY)^/^. 
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The saxion decay takes place far before the beginning of the big-bang neucleosynthesis 
(BBN), and hence is cosmologically harmless. It should be noted that the saxion might 
decay dominately into two axions. If it is the case, the extra energy of the axion at the 
BBN epoch raises the Hubble expansion of the universe, which leads to overproduction of 
^He. To avoid this difficulty the branching ratio of the saxion decay into two axion should 
be smaller than about 0.1. Here, we simply assume that the dominant decay process of 
saxion is a "^g-^ Then the saxion decay increases the entropy of the universe by the 
rate 

Because of this entropy production by the saxion decay, the primordial lepton asymmetry 
shown in Eq. (^Sj) is also diluted by the rate A. Then the present baryon asymmetry is 
given by 



4 X 10-^ <5eff- 



1/2 




Notice that the present baryon asymmetry is independent on the reheating temperature 
Tfj. We see that the desired baryon asymmetry is just given by the lightest neutrino 
mass of m^i ~ 10"^ eV for F„ ~ 10^° GeV. 

Before closing this section, we should comment on the cosmological consequence of 
our model. The entropy production by the saxion in Eq. (^) ensures that we are free 
from the cosmological gravitino problem. The number density of gravitinos produced at 
the reheating process is diluted by the rate A and hence the radiative decay of gravitino 
dose not disturb the BBN. Furthermore, we find that the model does not suffer from 
the cosmological problem of the axino which is a fermionic superpartner of the axion. 
The axinos are produced in the reheating process by the thermal scatterings and may 
lead to a cosmological difficulty p2[. However, in our model, the interaction of axino at 
the reheating epoch is suppressed by Xq ~ M^,, not by F^, and hence the production of 
axino is less effective. Furthermore, the entropy production by the saxion decay dilutes 
the axino abundance and hence the axino becomes completely cosmologically harmless. 



This is realized when m\ ~ m— . 
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6 Conclusions and discussion 



In this paper we have performed a detailed analysis on the Affleck-Dine leptogenesis 
taking into account the thermal effects from the dilute plasma. We have fflst shown that 
the thermal effects change drastically the dynamics of the flat direction and suppress 
the lepton-number asymmetry produced by the flat direction cj). In order to escape from 
the thermal effects, a relatively low reheating temperature of Tr ^ 10^-10^ GeV is found 
to be required when the lightest neutrino mass is niyi ^ 10~^ eV. This is consistent with 
the result obtained in Ref . . On the other hand, we have found that this upper bound 
on Tr is relaxed as Tr ^ 10^-10'^ GeV x(mjy]^/10 ^ eV) ^ for the lighter neutrino mass 
region of m^i ^ 10^^ eV. We have also estimated the resultant lepton asymmetry for the 
higher reheating temperature region where the thermal effects are important by both 
analytical and numerical calculations. We have found that an ultralight neutrino with 
a mass such as mj^^ ~ 10~^ eV is required to produce enough lepton asymmetry to 
account for the baryon asymmetry in the present universe. Here, we have assumed that 
the reheating temperature should be lower than about 10* GeV to avoid the cosmological 
problem of the gravitino of mass — 100 GeV-1 TeV. Such an ultralight neutrino 
seems to be very unlikely the case, since the recent Superkamiokande experiments @, 12 



suggest the masses of heavier two neutrinos z/2 and z/3 to be in a range of 10~^-10~^ eV. 

However, in the second part of this paper, we have pointed out that the above neutrino 
mass TTLyi ~ 10^* eV is not necessarily the mass to be observed today, if the heavy 
Majorana masses of the right-handed neutrinos are dynamical valuables in the early 
universe. We construct a model based on the Peccei-Quinn symmetry to demonstrate 
our point. We find in this model that the neutrino mass in the true vacuum can 
be as large as niyi ~ 10^^-10"^ eV to obtain ub/s ^ lO^^'^-lO^^^ for the reheating 
temperature of Tr ^ 10® GeV. Therefore, we have shown that the AD leptogenesis works 
well as long as the masses for the right-handed neutrinos are dynamical valuables in the 
early universe. 
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